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ABSTRACT 



The Transonic Turbine Test Rig installed in the Turbo-Propulsion 
Laboratory at the Naval Postgraduate School was designed such that 
losses in nozzles and rotor blading and overall performance charac- 
teristics can be determined for operation at several combinations of 
blade configurations, pressure ratios and speeds. This study relates 
the empirical data for stator loss coefficients to those determined 
for the converging-nozzle stator presently installed, and discusses 
reasons for some of the discrepancies. Also discussed are results of 
the flow nozzle calibration, locked rotor data, and temperature and 
inlet swirl effects. A method is proposed to reduce these effects 



for future work on the Turbine Test Rig. 
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2 
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T Turbine 

th Throat 

u Peripheral or tangential direction 

0 Stator entrance properties 

1 Stator discharge properties 

2 Rotor discharge properties 



11 



I. INTRODUCTION 



Much of the design and analysis of modern turbines is based on 
the experience of the designer and his use of empirical data available 
in the literature. As a re-»>3t, many investigations of test turbines 
are being carried out in an effort to more fully explain and mathe- 
matically model the flow phenomena in stationary and rotating channels. 
Present methods for design of axial flow machines base calculations 
on a reference diameter, usually taken as the mean of the hub and tip. 
A major simplification is then made by assuming that at any point in 
a cross-section, such parameters as total pressure and temperature and 
axial velocity have the same values as at the reference diameter. The 
overall characteristics of a stage are then determined, but depend on 
how accurately flow discharge angles and loss coefficients are known. 
Thus, it is extremely important to be able to predict these values as 
functions of blade geometry, Mach number and certain other flow para- 
meters. Such a method as just described is conventionally known as 
a ’’One-Dimensional Analysis”. 

The Transonic Turbine Test Rig (TTTR) located at the Turbopropul- 
sion Laboratory is a particularly desirable test machine because it 
permits evaluation of performance at many combinations of operating 
conditions and configurations. One of the most important features 
of the TTTR is that it enables one to evaluate flow angles and loss 
coefficients using momentum and moment of momentum equations, without 
making extensive flow surveys which cause an upset of the true flow 
conditions. The empirical methods available for prediction of flow 
angles and loss coefficients yield results varying from unreasonably 
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low values to one so conservative it would be unlikely the machine 
would operate. Results from previous runs on the TTTR showed the 
latter trend in loss coefficients, although the TTTR had a reasonable 
overall performance. 

This thesis investigates the methods by which the loss coefficients 
are determined for the TTTR and relates these coefficients to those 
predicted by the literature. It then describes the investigation car- 
ried out to explain the apparent discrepancies, and examines the effect 
of some of the more important parameters used to evaluate the loss co- 
efficient . 

The author would like to express his appreciation for the compe- 
tent assistance provided by Mr. Jim Hammer and members of the Turbo- 
propulsion Laboratory. Special thanks are extended to Dr. M. H. Vavra 
for his guidance and advice during the investigation. 
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II. EQUIPMENT MODIFICATIONS 



The TTTR instrumentation has been described extensively by Commons 
[Ref. 4] and Lenzini [Ref, 6). Discussion here will be limited to 
recent modifications made to the TTTR. 

Fig. 1 shows the piping and valve location in the Turbine test 
cell. Fig. 2 shows a cross-section of the floating stator and sup- 
port assembly, while Fig. 3 shows a more detailed cross-section of 
the stator plenum and hub area. The closure plate assembly shown in 
Fig. 3 was instrumented with additional strain gauges on the flex- 
ures, and the back face towards the rotor had three static pressure 
taps installed at known radii. Results of these pressure readings 
and closure plate torque calibration are discussed later. 

Fig. 4 shows the floating stator assembly with the locations of 
the axial force and torque capsule supports and Fig. 5 shows a more 
detailed view of the torque capsule support assembly. The capsule 
support assembly was moved from the previous location between the 
overhead support and plenum casing to its position at the upstream 
end of the stator in an attempt to reduce both the temperature effect 
on the capsule and the errors induced by differential thermal expan- 
sion. 

Fig. 6 shows £he stator exit and shroud area. Two additional 
pressure taps, P13 and P14 were installed, with the hope of corre- 
lating the stator discharge tip pressure to the shroud pressure at 
the same plane. In addition, the machine was completely dismantled 
and cleaned, and new pressure lines were installed in many cases. 

The pressure lines were connected to a mercury manometer board in a 
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logical sequence, extending from the nozzle to the plenum assembly 
and along the shroud to the rotor discharge, so that the pressure 
distribution through the machine could be readily observed. 

Fig. 7 shows the converging nozzles presently installed in the 
TTTR and Fig, 8 shows the circular-arc, sharp leading-edge rotor used. 
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III. FLOW RATE DETERMINATION 



A. NOZZLE CALIBRATION 

Previous nozzle calibration techniques have been described by 
Naviaux ^ and Lenzini [Ref. 6j. It was shown that the nozzle coeffi- 
cient was a function of nozzle Reynolds number. However, former re- 
sults for nozzle coefficients were determined as a function of the 
Reynolds number based on the conditions of the ASME standard calibrat- 
ing orifice, rather than on the flow nozzle. As a result, it was de- 
cided to repeat the nozzle calibration to obtain a new relationship, 
and to determine if any deterioration of the nozzle had taken place. 

The program for the data reduction is given in Appendix D, program 

2 

FLOCAL. The method used is similar to that given by Williams , which 

was found to be more suitable for computer computations. The basic 

3 

equation is given in the ASME Power Test Code and values of the con- 

4 

stants are taken from Stearns. Tests were carried out at nozzle sup- 
ply pressures of 26, 29, 34 and 39 psia over a range of Reynolds 
numbers from 2 x 10^ to 12 x 10^. Fig. 11 shows the results of these 



Naviaux, J. C., '’Transonic Turbine Test Rig Exhauster System Tests 
and Tests of a Reaction Turbine, "Naval Postgraduate School Thesis, 
December, 1966. 

2 

Williams, D. D., "Determination of Performance Parameters of a 
Dual-Discharge Radial Turbine, "Naval Postgraduate School Thesis, 
December, 1968. 

3 

Flow Measurement , Chapt. 4, Part 5, Supplement to ASME Power Test 
Codes, ASME, New York, N. Y., 1959, p. 57ff. 

4 

Stearns, R. F., et_ al. , Flow Measurement With Orifice Meters . 

New York: D. van Norstrand Company, 1951. pg. 65-212. 
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tests, with nozzle coefficient plotted as a function of Reynolds 
number. The curve differs from that given by Lenzini [Ref. 6], being 
higher by 2 to 3 percent. The curve also shows a maximum coefficient 
occuring at a Reynolds number of 6 x 10^ and then decreasing as Rey- 
nolds number is increased. It is suspected that this increasing trend 
of nozzle coefficient each time a calibration is performed indicates 
some deterioration of the nozzle. 

An analytical curve for the nozzle coefficient was obtained by us- 
ing a least-squares polynomial curve fit. A fourth-order polynomial 
was chosen as the best fit to the data, namely: 

C = 9.5208828 x 10 _1 + 4.4891640 x 10” 7 Re -7.5810743 x 10” 13 Re 2 
n 

+ 5.5023778 x 10” 19 Re 3 -1.4571314 x 10" 25 Re 4 

The maximum error between the analytical curve and the actual data 
points for a range of Reynolds numbers between 4 x 10^ and 12 x 10^ is 
0.2 percent. 

B. TURBINE FLOW RATE 

The program to obtain the turbine flow rate is given in subroutine 
FLORAT of program TTTR, and employs a procedure similar to that de- 
scribed in section A. The labyrinth leak rate uses a formula deter- 
mined by Lenzini [Ref. 6] which is a function of plenum pressure ratio 
and temperature. The labyrinth leak rate was found to be only 2 per- 
cent of the total flow rate. The recalibration of the nozzle and im- 
proved flow rate determination resulted in referred turbine flow rates 
which were 3 percent higher than those previously used. 
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IV. ANALYSIS AND DATA REDUCTION 



A. GENERAL 

The basic equations for a one-dimensional flow analysis of the 
TTTR are well known and have been described in detail by Commons 
[Ref. 4], Lenzini [Ref. 6] and Messegee [Ref. 7]. The reference diam- 
eter referred to previously is twice the mean radius: 

\ub + R tip 

R — o 

m 2 

The data for the TTTR are reduced by program TTTR given in Appendix D, 
which also gives the results of Run 15. Fig. 12 shows the thermody- 
namic process for an axial turbomachine and the stator loss coefficient 
is designated Fig. 13 shows the associated velocity diagram for 

a turbine stage. 

It is of interest here to summarize the assumptions made in using 
the one-dimensional analysis, as they will have a large bearing on 
the accuracy of the results. These assumptions are: 

1. Adiabatic, steady flow, so that the absolute and relative en- 
thalpies are constant respectively along any streamline. 

2. Cylindrical stream surfaces, so there is no bending of the 
stream surface, or radial component of velocity. 

3. Uniform conditions ahead of the stator so there is no total 
pressure or total temperature gradient in the radial direction. 

4. Completely radial flow into the plenum, and axial flow into 
the stator with no tangential velocity component. 

5. Shear stresses are negligible on the annulus and shroud, and 
symmetrical wakes occur at the trailing edges of the blades. 
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The validity of these assumptions will be discussed later and 
where invalid assumptions were found to exist their . inf luence on the 



loss coefficient will be discussed. 

When analysing the reduced data, nearly all variables were plotted 
as functions of the stator discharge pressure ratio P^/P^. The 
following discussion will show why an i m P ort:ant: parameter. 

From the equation of momentum 

w 

M = 77 — RV 
g c u 



and using 



= 



w T vV 8 



g & c TO 



P T0 A ths 




1 -I 



n - 1 
n 



TO 



where 



n = 



1 + c p (Y - 1) 



and 



C_ = 



P - P 
TO T1 

P - P 
TO 1 



hence 



4> = f 



TO 



1 l P, 



, n, y , 



Thus 



M = 



P^ A , 4> 

TO ths 



n/V* 



T 

c TO 



RV 



or 



JL . H . \ 

ref /=— 



TO 



TO 



TO 

'll P 1 * 



n, Yj 
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Now 



„ 2 2 . 2 
V = V sin a 
u 



- C ( T I0- T l) sl 



2 

sin a 



or 



V U =CT X0 



1 - 



(V\^ 

VIC'/ _ 



. 2 

sm a 



and 



\A 



= V* 



Hence 



TO 



M = f 



1 - 



n - 1 
n 



TO 



sin a = f 



TO 



2 \ P„ 



n, Y/ 



TO 



ref 3 V P, 



n > Y * a 



For an ideal gas y is constant, and for small pressure ratios n is 
nearly a constant. Thus, if it is assumed that the stator discharge 
angle is constant: 



M ref f 4 { P 



TO 



and 



and 



$ = f 



TO 



1 V P. 



V 



= f 



TO 



\/ t to 



2 P, 



Using the loss coefficient as defined by Vavra [Ref. 8]: 



■ 1 - » 2 - 1 - ( V 1 



V, \ 2 T - T 



T /T - T /T 
lis 1 TO lis 7 TO 



lis 



T - T 
TO lis 



1 - T /T 

lis 7 TO 



n - 1 



( P 1 /P tJ - (y P T^ 

? ]/ P To) 



Y ~ 1 
Y 



Y - 1 



= f 



TO 



s l P, 



n> Y 
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And, consistent with the above assumptions for n and , 




B. DESCRIPTION OF METHOD 1 



It is possible to solve for the stator discharge properties by 
satisfying continuity at the throat plane and the exit plane, provided 
certain parameters are accurately known and some assumptions are made. 
At the exit plane, from continuity: 



“T * C 1 V al A 1 



From the perfect gas relationship 

P, 



P 1 R T, 
g 1 



Thus 



al 



W T R g T 1 
p i A i 



(i) 



However T^ is not known until is determined. If one assumes that 
the same pressure and temperature acting at the exit also exist across 
the throat plane, then continuity may be applied at the throat: 



“T ■ P 1 V 1 A ths 



Then using 



V. 



p_ = - — — and T_ = T__ 
M 1 R T„ 1 TO 



8 . 1 
w_ 



V, = 



“ T T 1 



1 P, A . P, A . 

1 ths 1 ths 



w T R 

-I . & 

p i A , 

1 ths 



or 



o C P, A , 

V 2 + — — V - 

1 - „ 1 P, A 



“t r e 



"t R g 



1 ths 



TO 
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Let 



C 



2 



C P, A , 

1 ths 



w n 



R 

g 



Thus 




C T 



TO 



And completing the square: 




( 2 ) 



(3) 



Once V i is determined from (1) , 
al 



= cos 



-1/ V al 



= cos 



-1 



ths 1 



and V , = V- sin a_ 
ul 1 1 



1 / V*l 

and all quantities are available to determine £ . 

s 

It is important to examine in detail the assumptions made in 
arriving at by continuity only, which are: 

1. The same pressure P^ and temperature T^ exist at the throat 
and exit plane, 

2. There is no change of from the throat to the exit plane. 

3. The correct pressure distribution for is accurately known. 

4. The areas A . and A- are known exactly. 

ths 1 J 

5. There is one-dimensional or uniform flow at both the throat 
and exit plane, which implies no boundary layer restriction. 



In the data reduction, P^ was taken as the mean of the measured 
hub and tip pressures, and A^ was the annulus exit area reduced by a 
factor K to account for the trailing edge thickness. 
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C. DESCRIPTION OF METHOD 2. 



In this method both continuity and momentum are applied, where 
the tangential velocity is determined from the momentum equation, 
while the axial velocity is determined from knowledge of the axial 
forces acting on the stator and by applying continuity at the stator 
exit. From angular momentum: 

M g 



ul 



”t r „i 



( 4 ) 



At the stator exit plane, a parabolic pressure distribution is assumed 
to exist, which is expressed by equation 12 of Lenzini [Ref. 6] and is 
repeated here: 

P, 



hub 



(1 + R tip + R tip “hub - < 2 + "hub 



lav 



(R 



tip 






P . 

+ -^ 



(R tip - "hub 2 ) 



(5) 



The value of e determines the shape of the parabola. A value for e 
is first assumed and with this first approximation for P^^, the force 



balance equation may be solved for V , 



V 



alM 





( 6 ) 



where R = the algebraic sum of the measured forces and the pressure 
forces acting on the stator. Thus V ivr is a direct function of the 
value of e, or the shape of the pressure distribution. Continuity is 
then written for the exit plane as in Method 1, employing the perfect 
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gas relations, to obtain: 

v .vi ( T .!i! 

alC P A 2 \ TO C 



W T \ 

P 1 A 1 



TO 



1 

C 



V ^ + V 2 

alC ul 



Let 



c = W T R g 
1 P 1 A 1 



Then 



C 2 C 

V alC = C 1 T T0 - - V alC ' “ V ul 



or 



V alc + C L V alc CT T0 



and completing the square: 



V 



2 

ul 




or 



V 



ale 




Equations (6) and (7) are compared and the value of e is changed appro- 
priately until the two methods produce the same value of V In sum- 

mary then, the steps are: 

1. Assume an initial value of e and determine P from equation 

(5). 

2. Solve for from equation (4). 

3. Solve for V ^ from equation (6). 

4. Solve for from equation (7). 

5. Compare the two values of V ^ and vary e until the two values 

of become equal. 

6. Compute £ . 



25 



The assumptions made using this method are as follows: 

1. The stator torque can be measured accurately, yielding a cor- 
rect value of V - . 

ul 

2. The correct pressure distribution for is known accurately, 
or conversely the iterated value of e will yield the correct pressure 
distribution. 

3. The correct axial force is being recorded. 

4. The area is known exactly. 

5. There is one-dimensional or uniform flow in the nozzles, so, 
no boundary layer correction factor is required. 

As in Method 1, the area A^ was reduced by the same factor K._ to 

account for the effect of trailing edge thickness. 

If the assumptions for both methods are met, then either method 

should result in essentially the same values for V - , V _ , a, and L , 

al ul s 

respectively . 
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. V. NATURE OF PROBLEM 



A survey of previous results has shown that C g has ranged from 
negative values to as high as 0.30. The result of this has produced 
stator or rotor efficiencies of over one-hundred percent, but the 
overall machine had a total-static efficiency in the order of eighty 
percent at the design point. The conclusion reached was that the 
determination of overall performance was correct, but that some error 
existed, either in the theory of calculating stator loss coefficients 
or in measuring the properties necessary to arrive at the loss coef- 
ficient . 

Appendix A presents several methods by which the stator discharge 
angle and stator loss coefficients may be predicted. In addition, a 
curve fit of data is given by Vavra [Ref. 8]. A summary of these 
values and some representative ones from previous work by Lenzini 
[Ref. 6] are presented in Table I below. 



TABLE I 

Predicted and Experimental Values of £ 

s 





Source 


c 






s 




Ainley 


0.0756 


Predicted 


Traupel 


0.1545 


Values 


Soderberg 


0.1022 




Vavra 


0.1400 


Experimental 


Lenzini 

j 


-0.0100 to 0.3000 

| 
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Since previous work by Lenzini [Ref. 6] indicated that the over- 
all total-static turbine efficiency was about eighty percent, it was 
concluded that many of the above values were too high, and that an 
average value for the loss coefficient of 0.08 would be more in 
order. The problem then became one of determining where the error 
in computing the loss coefficient was occuring and to reduce the error 
as much as possible. 

In order to determine the error, it was assumed that the theory 

to calculate C was correct and that the error arose from invalid 
s 

assumptions and incorrect measurements of some variables. If a value 
for is assumed, certain parameters can be calculated and plotted 
as a function of and then actual measurements could be com- 

pared against these predicted values. Thus, based on the definition 
of > several important relationships could be developed: 




Let X - 1 



Hence 



and using C = 





- X 




/ CT 



) 



i ■ x/Ht 0 - OC 1 - *) 

With the chosen value of z; = 0.08 and y = 1.4, 



( 8 ) 
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— = 2.1450 (1 - X) 2 
a T0 
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VI. DESCRIPTION OF TESTS 



A new electrical read-out system was employed for the strain- 
gauges on the closure plate. It was decided to calibrate both the 
closure-plate axial force and torque with the closure-plate mounted 
on the stator, rather than using the previous method of calibrating 
away from the TTTR and then possibly destroying the calibration during 
the installation in the stator casing. 

A summary of the operating runs is listed in Table II. The first 
series of runs was carried out at pressure ratios of 1.5 and 2.0 over 
the full range of RPM allowable by the torque capacity of the air dy- 
namometer. The configuration was unhooded with rotor axial clearances 
of 0.250 in. and 0.500 in. Several pressure leaks developed on nearly 
every run, and only Run 8 was considered sufficiently trouble-free to 
be used for analysis. 

Once the high loss coefficients were confirmed to exist from the 
data of Run 8 and reduction by Method 1 and Method 2, it was decided 
to attempt some locked rotor tests at axial clearances of 0.250 in. 
and 0.500 in, and 2.000 in. to determine the conditions at zero RPM 
and a pressure ratio of 1.5. Few problems were incurred during this 
series of tests, and data points could be repeated to within 1 percent 
or better. During the last run however, several pressure lines were 
torn loose, which required extensive repairs before any further runs 
could be made. 

During the runs up to and including Run 14, the standard method 
of operation had been to leave the inlet valve of Tank 1 full open and 
gradually open the turbine inlet valve to increase the pressure ratio. 
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TABLE II 



Summary of Runs 



Run 


Configuration 


Purpose 


1-4 


ASME Standard Pipe and Orifice Installed 


To determine nozzle 




P varied 


calibration coefficient. 




noz 




5 


Rotor removed 


Closure plate axial and 


6 


Rotor removed 


torque calibration 


7 


Operating, PR = 2.0, AX = 0.250" 


To check the system 


8 


Operating, PR = 1.5, AX = 0.250" 


after modifications 


9 


Operating, PR = 1.5, AX = 0.500" 


and determine operating 


10 


Operating, PR = 2.0, AX = 0.500" 


data 


11 


Operating, PR = 2.0, AX = 0.500" 




12 


Locked Rotor, AX = 0.500" 


To examine the system 


13 


Locked Rotor, AX = 2.000" 


operation at zero RPM 


14 


Locked Rotor, AX = 0.250" 


and various clearances. 


15 

| 


Operating, PR = 1.5, AX = 0.250" 


To confirm previous 
results 


16 


Honeycomb, flexures installed 


To determine the ef- 


17 


Honeycomb, flexures removed 


fects of the torque 


18 


Honeycomb, flexures installed 


flexures and various 


19 


Honeycomb , flexures removed 


turbine inlet valve 
openings . 


20 


Nozzles installed, rotor removed, 


To determine the ef- 




flexures removed 


fects of torque flex- 


21 


Nozzles installed, rotor removed, 


ures and inlet valve 


J 

f 


flexures installed 

! 

I i 

1 ! 

j 


openings with no rotor 
(influence. 
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As a result, a pressure drop in the order of 50 in. of mercury was 
occuring across the turbine inlet valve. Also, it was noticed that 
after shutting off the flow to the stator, there were residual values 
of stator axial force and torque in the order of 1.0 and 5.0 percent 
respectively, which gradually dropped to zero after some time. 

In an effort to examine the validity of the assumption of no in- 
let tangential velocity ahead of the stator, the stator was removed 
and the honeycomb flow straightener shown in Fig; 9 was installed. 
Also, during the next series of tests, the method of operation was 
changed and a further modification was employed. '-For Run 16, the 
honeycomb was installed and the turbine valve operated as before. 
Because a residual torque reading was again observed, it was decided 
to remove the stator torque flexures which take up some of the torque 
and transmit the axial force to the cradle assembly. A system of 
roller bearings was installed between the cradle assembly and stator 
plenum casing in place of the torque flexures, as shown in Fig. 10. 
Also, a thrust collar was fastened on the stator support shaft and 
allowed to press against the support bearing, as a means of transmit- 
ting the stator axial force to the cradle assembly. The next series 
of runs was performed by either gradually opening the turbine valve 
to obtain the desired pressure ratio, or fully opening the valve and 
adjusting the pressure ratio by slowly closing a dump valve in the 
delivery line near the Allis-Chalmers compressor. It was possible 
to obtain similar flow rates as used during all previous tests, but 
because of the much less restricted flow, considerably smaller pres- 
sure ratios were required. 
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Prior to Run 16, a thermocouple was installed on the torque 
capsule support assembly, and during the next series of tests the 
temperature of the assembly was monitored. Also, the discharge tem- 
perature from the Allis-Chalmers compressor was varied from 90° to 

¥ 

130° F to provide a sufficient temperature gradient at the stator in 
order to examine the effects o.i differential thermal expansion. Dur- 
ing Run 17, the flow exit angle from the honeycomb was measured using 
a 5-hole probe, to find out whether the flow leaving the honeycomb 
straightener was truly axial as desired. 

Once an improved method of operation was determined, a final 
series of tests was performed with stator nozzles reinstalled, but 
the rotor removed to eliminate any influence it might have. The same 
procedure was adopted as above, with the flexures installed and the 
turbine valve operated either fully open or partly open, and then the 
flexures were removed and the roller bearing supports installed, and 
the valve operated either fully open or partly open. 
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VII. RESULTS AND ANALYSIS OF TESTS 



Fig. 14 shows the total-static turbine efficiency as a function 

of isentropic head coefficient. It is seen that a maximum efficiency 

of 79.6 percent occurs at a k. of 4.2. Previous results for these 

is 

operating' conditions quoted a maximum efficiency of 83 percent at the 

same k. . However, the efficiency is directly proportional to the dy- 
1 s 

namometeY torque M, and inversely proportional to the flow rate w^ 

d 1 

for given inlet conditions and speed N. Since the calculated flow 
rate is now considered to be in the order of 3.0 percent higher be^ 
cause of the corrected nozzle coefficient, this would account for the* 
difference, and it was therefore concluded that the correct results 
were being obtained. 

Fig. 15 shows the referred stator axial force as a function or 
referred RPM. It can be seen that there is some non-linearity when 
the curve is continued to the locked rotor condition, which can be 
accounted for by the changing pressure distribution around the shroud 
with the rotor locked. 

Fig. 16 shows the referred stator torque as a function of referred 
speed. There is good linearity to the locked rotor condition, indi- 
cating the locked rotor torque is in the range of its expected value. 
When the referred torque is plotted as a function of stator pressure 
ratio, as shown in Fig. 17, there is a difference between the results 
of Method 1 and Method 2 and a considerable difference between the 

actual and the predicted value for the assumed £ =0.08. The dif- 

s 

ference between the methods is due to the resulting higher pressure 
of P^ by Method 2. It can be concluded however, that the torque 
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measurements at a pressure ratio of 1.45 are about 8 percent too 
high. It will be discussed later how this can occur. 

Fig. 18 shows the flow function $ plotted against the stator pres- 
sure ratio. The curve extends well to. the locked rotor condition. 

The reference curve shows that the actual $ is more in error at the 
lower than at the higher pressure ratios when using Method 1. The 
difference between the curves of Method 1 and Method 2, as discussed 
previously, is due to the higher by Method 2. Although $ is a 

function of both w^ and A- for fixed inlet conditions, it was con- 

T ths 

eluded that A^g was in error. Build-up of dirt deposits on the 
blades and the presence of an increasingly thick boundary layer at 
higher pressure ratios would result in a lower effective throat area 
than was used. 

Fig. 19 shows the dimensionless velocity ratios V^^/a^ and 

^ul^ a T0 aS ^ unct ^ ons the stator pressure ratio. In both cases 

the curves extend smoothly to the locked rotor condition, and there 

is generally good agreement between the two data reduction methods. 

It can be seen that both V - /a mr . and V - /a^ are lower than the cal- 

al 1U ul J.U 

culated values for an assumed £ =0.08. When using Method 1, V - 

s ax 

is determined from equation (1) , and a low value of V . could be due 

al 

to a high value of pressure P^ or the reduction factor K^. When us- 
ing Method 2, V ^ is determined by varying e until equations (6) and 
(7) agree, which means in effect that P^ is varied. Method 2 implic- 
itly assumes that V is correct. Examination of the behaviour of 
equations (6) and (7) when reducing the data revealed that equation 

(7) generally produced a lower V than equation (6), and its value 

al 

changed very little during the iteration. The strongest variables 
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in equation (7) are and K^_ , and a low value of would be due 

to high values of P and , which is in agreement with the previous 

conclusion for Method 1. The magnitude of V . has little effect on 

V ^ in equation (7) . Since by equation (6) is generally higher 

than by equation (7), and because changing P^^ causes equation (6) 

to change more rapidly than equation (7), the iteration scheme forces 

P^ higher in order to reduce the net stator force and thus reduce 

in equation (6) . The increase in P^ causes only a slight reduction 

in V 1 by equation (7). If the initial value for P- with e = 0 were 
al i 

too high then the above undesirable situation would occur, and the 

net result would be too high a value of P- and too low a value of V . 

1 al 

A desirable condition would occur if the initial value of by 

equation (7) were higher than that by equation (6), since the iter- 
ation scheme would then reduce P^ causing a higher resultant V^. 
Previous work by Commons [Ref. 4] indicated that the stator hub pres- 
sure and the stator tip pressure measured with the arrangement of 
Fig. 6 were 5 percent higher and 5 percent lower, respectively, than 
the pressured obtained by flow surveys after the stator. The over- 
all effect of these errors would produce a higher initial P^ than the 
correct value, which is consistent with the previous conclusion for 
the low value of V^. Examination of the results from Run 15 given 
in Appendix D shows that the hub pressure acting on the hub area pro- 
duces the largest effect on the net force in equation (6) . As a con- 
sequence, a slightly low value of P would result in the undesirable 
condition of a higher value of V ^ by equation (6) than by equation 
' (7). 
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When using Method 1, the low value of v ul / a T Q i n relation to the 

calculated curve for an assumed L =0.08 would be due to the low 

s 

values of and caused by an erroneously high pressure P^. This 

supports the previous argument where it was stated that high values 
of P^ could be causing the errors. Using Method 2, equation (4) shows 
that the low value of V ^ could be caused by a low stator torque read- 
ing. However, that condition would conflict with the conclusion from 
Fig. 17 which shows the torque to be reading about 8.0 percent high. 

A reduction in the stator torque would result in a further deviation 

of V from the calculated curve for £ =0.08, and an increase in 

ul s 

the mean loss coefficient. This situation can be resolved only by 
concluding that the angle of 75.9° chosen to calculate the velocity 
components from equation (9) was too high. A reduction in the chosen 
flow angle would result in a higher and lower for a given V^. 

Fig. 20 shows the stator loss coefficients calculated by the two 
methods as a function of stator pressure ratio. It can be seen that 
the losses obtained by either method are at least 50 percent higher 
than the predicted value. The higher loss coefficient by Method 2 
is due to the higher pressures resulting from Method 2 and the fact 
that constant pressure lines diverge. Fig. 20 also shows that the 
loss coefficients determined by both methods decrease slightly as the 
stator pressure ratio and the flow rate increases. Such a condition 
is opposite to that expected, since at higher pressure ratios an in- 
creased possibility of separation exists, and a thicker boundary layer 
with greater secondary losses should occur. The apparent reduction of 
the loss coefficient as the stator pressure ratio increases is likely 
due to an unexpected flow condition at higher flow rates. It will be 
shown later that at higher flow rates a changing inlet flow condition 
does occur. 
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Fig. 21 shows the referred stator axial force as a function of 
stator pressure ratio for three axial clearances of 0.250 in. 0.500 
in. and 2.000 in. It can be seen that at a clearance of 0.500 in. a 
more negative axial force occurs than at 0.250 in. or 2.000 in. It 
is suspected that the presence of the locked rotor at various axial 
clearances causes a changing pressure distribution on the shroud, and 
thus the rotor axial clearance has a large effect on the stator axial 
force and the shroud pressure distribution. 

Fig. 22 shows the referred stator torque as a function of the 
stator pressure ratio. The same reversal of trend occurs when the 
axial clearance is increased, confirming the fact that the presence 
of the rotor influences conditions upstream. The locked rotor results 
shown in Fig. 22 at the three different axial clearances further sup- 
ports the previous conclusion from Fig. 17, namely that the stator 
torque is reading high. 

Fig. 23 shows the referred flow rate as a function of stator pres- 
sure ratio for the three axial clearances at the locked rotor config- 
uration. It can be seen that there is good agreement at each of the 
three axial clearances, so the changing torque with clearance is not 
due to a different flow rate. Fig. 23 also shows that the flow rate 
approaches a maximum as the pressure ratio nears the critical value, 
but that choking does not occur. 

Fig. 24 shows the flow function $ with changing stator pressure 
ratio for the locked rotor condition. As in Fig. 23 there is no ef- 
fect of clearance on flow rate. However, the same trend of a lower 
flow function than desirable is occuring, which confirms the results 
of Fig. 18. 
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Fig. 25 shows the stator discharge pressure ratio as a function 
of the overall pressure ratio p T q/ P 2 * It: can be seen there is a 
linear relationship between stator pressure ratio and overall pres- 
sure ratio. Fig. 25 also shows the increasing difference between 
Method 1 and Method 2 as the overall pressure ratio is increased. 

Fig. 26 shows the referred velocities as a function of stator 
pressure ratio for the locked rotor configuration. The same trend 
in V ^ occurs with the locked rotor results as was noticed in Fig. 19 
when the rotor was unlocked, confirming the conclusion that some er- 
ror exists in the pressure and the reduction factor . However, 
Method 2 results in a tangential velocity component which is greater 

than the calculated value for the assumed £ =0.08. It can there- 

s 

fore be concluded that the increased value of V ^ is directly 

ul TO 

due to an error in the stator torque measurement and that the high 
torque reading is caused by some combination of temperature effect 
on the assembly and an inlet swirl at the stator plenum. 

Fig. 27 shows the stator loss coefficient as a function of stator 
pressure ratio for the locked rotor condition. Equation (8) gives 
the relationship between the pressure ratio, exit velocity and stator 
loss coefficient, and it can be seen that for a given pressure ratio 
a low velocity will result in a high loss coefficient, and conversely, 
if the velocity were too high, the loss coefficient could become neg- 
ative. When using Method 1 the loss coefficient is a unique function 
of the throat pressure and area where the throat pressure was assumed 
to be the same as that existing at the exit plane. Behaviour of equa- 
tion (2) shows that a low velocity will result if a combination of 
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the pressure and throat area is too large. It has been previously 
shown that the throat area appears to be too large and Figures 20 and 
27 show that at low pressure ratios the throat pressure must be high 
also. Appendix C gives a predicted pressure distribution at the exit 
plane, and it can be seen that a parabolic, rather than a linear var- 
iation would be more in order, resulting in a slightly higher than 

the linear variation from hub to tip would predict. From this it may 
be concluded that assumption 1 of Section 4.B, namely that equal pres- 
sures exist at the throat and exit plane would be in error. As the 
stator discharge pressure ratio increases the loss coefficient de- 
creases somewhat, indicating that the pressure at the throat and exit 
plane are more nearly equal. When using Method 2, an erroneously high 
stator torque reading will produce a high tangential velocity compo- 
nent resulting in a higher net velocity than should exist for a given 
pressure ratio. It should be noticed that the loss coefficient by 
Method 2 increases with increasing pressure ratio even though the 
stator torque is high. This indicates that the error in hub pressure 

P, , and exit pressure P- have a stronger effect than does the stator 
hub 1 

moment and it was concluded earlier from the analysis of Fig. 19 that 
P^ was indeed too high. 

Fig. 28 shows the measured temperature of the torque capsule sup- 
port assembly as the flow rate is increased, and Fig. 29 shows the 
results of the temperature analysis on the cradle assembly. It can 
be seen that as the cradle assembly heats up, a positive error occurs 
in the stator torque reading, indicating that the reading becomes too 
high as the inlet temperature increases. The reason for the increasing 
torque error is the unequal expansion of the capsule support assembly 
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and the connectors between the point of action of the torque arm and 
the upper end of capsule support. The stator torque flexures pre- 
vent free rotation of the stator as the capsule assembly grows due to 
thermal expansion. The two steel rods and the Invar ring of the torque 
capsule do not expand by an equal amount in the opposite direction to 
cancel the thermal expansion of the capsule support. The differential 
thermal expansion causes an additional tensile force on the capsule 
and a higher apparent torque reading. Additional analysis showed that 
it would be impractical to compensate for this thermal expansion by 
making the support assembly of the same type of material as the capsule 
ring and steel attachments, since the compensation would be possible 
only for a fixed temperature. In addition it would be impractical to 
thermally shield or cool the assembly as was done with the dynamometer 
capsule, since this set-up would present severe complications if the 
turbine were operated inside the closed hood. Thus the only practi- 
cal way to eliminate the thermal effect on the capsule support assembly 
was to give the stator freedom to rotate so that no additional tension 
could be applied to the capsule. A solution was accomplished by re- 
moving the torque flexures and installing roller bearing supports to 
maintain the stator alignment, and a thrust collar was fastened to the 
stator support shaft and allowed to press against the support bearing 
as a means of transmitting the stator axial force to the cradle assem- 
bly . 

Fig. 30 shows the referred stator torques obtained with the honey- 
comb flow straightener replacing the guide vanes. Test data are shown 
at different openings of the turbine inlet valve and the tests were 
carried out with the torque flexures installed or removed. If the 



41 



flow after the honeycomb were axial, and if the flow entered the sta- 
tor assembly with no whirl components, the stator torque would be zero 
for any pressure ratio. Fig. 30 shows that the stator torque was only 
zero at zero flow rate, and increased with increasing pressure ratio. 
From the results of Fig. 30 it can be concluded that temperature ef- 
fects caused by the presence of the torque flexures and throttling the 
flow by partly closing the turbine inlet valve influence the stator 
torque readings. With the turbine valve full open less torque occurs 
when the flexures are removed, indicating that the temperature effect 
due to the flexures has been reduced. When the pressure ratio is set 
with the turbine valve partly closed a considerably larger torque 
occurs and it is again higher if the torque flexures are installed. 
Since an increasing torque with increasing pressure ratio still occurs 
when the turbine valve is full open and the torque flexures are re- 
moved, it can be concluded that either the honeycomb is not truly 
axial or the flow still enters the stator plenum with a small whirl 
component. 

Fig. 31 shows the flow angles measured after the honeycomb 
straightener . It can be seen that the velocity after the honeycomb 
has a tangential component which increases with increasing flow rate. 
From Figures 30 and 31 it can be concluded that the partly opened tur- 
bine valve causes whirl components in the delivery pipe, and the flow 
enters the stator assembly with a tangential velocity. Hence, assump- 
tion 4 of Section 4. a, namely, that the flow into the plenum is radial 
and no whirl component exists ahead of the stator is not satisfied 
with the present installation. Because the honeycomb has openings of 
about 1/8 inch it is apparently not capable of completely turning the 
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flow into the axial direction. It was found also that the stator 
blockage factor £ through the honeycomb is about 0.72, where £ rep- 
resents the ratio of the actual flow area and that necessary to pass 
the flow if the expansion were isentropic. Because of the low valves 
of 5 it can be concluded there must exist large whirl components ahead 
of the honeycomb causing flow separation at the inlet and large losses. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



From the preceding analysis, it can be concluded that the stator 
torque reading is about 8 percent too high, and that this error is 
due to a combination of the temperature effect on the capsule support 
assembly and an inlet swirl caused by the turbine inlet valve and 
piping. It is recommended that the turbine be operated with the tur- 
bine inlet valve full open at all times, and that the pressure ratio 
be set by closing a dump valve near the air supply. 

During the last series of runs with the stator nozzles installed 
and the rotor removed, inconsistent and residual values of stator 
axial force occured when the torque flexures were removed and the 
roller bearings installed. This was due to the mechanical friction 
inherent in the present system. It is felt that a better system 
should be designed which would allow the stator to rotate and move 
axially without restraints. This would make the stator a truly free- 
floating assembly and eliminate the influence of temperature on the 
measurement of the stator torque. Such a system would require the 
use of ball-bearings mounted in such a way that they would allow both 
rotation and axial movement of the stator and thermal expansion or 
contraction of the casing, but maintain the stator axial alignment. 

The analysis of the cradle assembly showed also that thermal ef- 
fects occur in the horizontal direction which affect the capsule meas- 
uring the stator axial force. This problem could be eliminated by 
constructing a cradle support assembly similar to that used for the 
torque measurement but made of aluminum. Proper installation of this 
cradle would allow its growth or contraction to be in the same 
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direction as the horizontal members of the stator cradle, thus com- 
pensating for thermal effects on these members. 

It is suggested that a honeycomb flow straightener and series of 
screens be installed in the delivery pipe just ahead of the inlet to 
the stator plenum chamber, and that the conical screen presently in- 
stalled in the plenum chamber be removed. This should produce a more 
uniform flow into the plenum and eliminate the tangential component, 
while removal of the conical screen should eliminate some of the 
existing total pressure gradient at the entrance to the stator noz- 
zles . 

From the analysis of the pressure effects on the loss coefficient 
it can be concluded that extensive flow surveys should be done at dif- 
ferent operating conditions to accurately determine the pressure dis- 
tribution at the stator discharge. Means must then be found to ensure 
that accurate hub and tip pressures can be recorded since they have a 
large effect on the loss coefficient. It was shown that the assump- 
tion of stator exit pressure also acting at the throat plane was er- 
roneous and the assumption could be eliminated if the throat hub and 
tip pressure taps were accurately located at the throat plane. At 
present there is up to a 5 percent error between the four sets of 
taps. The flow surveys must be carried out with probes which prop- 
erly compensate for Mach number effect, depth of immersion and pitch 
angle. Very useful information could be gained by using a hot wire 
anemometer probe to accurately determine velocity profiles and mag- 
nitudes in the stator area at different operating conditions. 

Results of the prediction of closure plate torque and shroud torque 
given in Appendix B show these values never exceed 0.1 percent of the 
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stator torque reading. This amount is considered negligible and 
thus the closure plate torque webs, strain gauges and wires could 
be removed. This would eliminate the congestion in the hub area 
and allow the closure plate greater freedom of movement, since some 
friction might be caused by the presently installed equipment. The 
predicted torque as a function of pressure ratio and speed agreed 
very well with the strain gauge readings and so it* is concluded the 
closure plate and shroud torques could be accounted for during the 
data reduction by using an empirical formula. Results of the forces 
on the closure plate due to the pressure distribution on the face 
agreed to within 1.0 percent of the strain gauge reading. Consid- 
erable effort was required to numerically integrate the pressure 
over the face area, and since measurement of the closure plate axial 
force by the strain gauges was much quicker and produced comparable 
results, it is recommended the closure plate pressure lines be re- 
moved . 

It can be concluded that several of the assumptions utilized to 
solve the equations for a one-dimensional analysis were not valid, 
and that ignoring their effects can produce very poor results. Appen- 
dix C shows that cylindrical stream surfaces cannot exist, and a 
radial component of velocity must be taken into account. Flow surveys 
made by Commons [Ref. 4] show that there is a variation of the total 
pressure at the stator entrance and thus a radial entropy gradient 
exists. The results of this investigation have shown that the flow 
is not completely radial into the plenum or axial into the stator 
and that a tangential component of velocity does exist. However, 
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these unsatisfactory conditions can be reduced by incorporating the 
equipment modifications suggested. 
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FIGURE 2 

TRANSONIC TURBINE TEST RIG 
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STATOR TORQUE CAPSULE SUPPORT ASSEMBLY 




FIGURE 6 

TURBINE AND SHROUD DETAILS 
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FIGURE 7 

MEAN RADIUS BLADE PROFILE CONVERGING NOZZLES 
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HONEYCOMB FLOW STRAIGHTENER 
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STATOR ROLLER BEARING SUPPORTS 
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FIGURE II 

NOZZLE COEFFICIENT VS. NOZZLE REYNOLDS NUMBER 




THERMODYNAMIC PROCESS FOR AN AXIAL 
TURBOMACHINE 
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VELOCITY DIAGRAM OF TURBINE STAGE 
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TOTAL- STATIC EFFICIENCY VS. ISENTROPIC HEAD COEFFICIENT 
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FIGURE 15 

REFERRED STATOR AXIAL FORCE VS. REFERRED RPM 
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FIGURE 16 

REFERRED STATOR TORQUE VS. REFERRED 
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FIGURE 17 

REFERRED STATOR TORQUE VS. STATOR PRESSURE RATIO 
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FIGURE 19 

REFERRED VELOCITIES VS. STATOR PRESSURE RATIO 
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FIGURE 20 

STATOR LOSS COEFFICIENT VS. STATOR PRESSURE RATIO 
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FIGURE 21 

REFERRED AXIAL FORCE VS. STATOR FRESSU RE RATIO 
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FIGURE 22 

REFERRED STATOR TORQUE VS. STATOR PRESSURE RATIO 
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FIGURE 23 

REFERRED FLOW RATE VS. STATOR PRESSURE RATIO 
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FIGURE 24 

FLOW FUNCTION VS. STATOR PRESSURE RATIO 
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FIGURE 25 

STATOR PRESSURE RATIO VS. OVERALL PRESSURE RATIO 
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REFERRED VELOCITIES VS. STATOR PRESSURE RATIO 
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FIGURE 27 

STATOR LOSS COEFFICIENT VS. STATOR PRESSURE RATIO 
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TEMPERATURE INCREASE OF CAPSULE SUPPORT VS. REFERRED FLOW RATE 




PREDICTED ERROR IN STATOR TORQUE DUE TO TEMPERA- 
TURE EFFECT ON CAPSULE SUPPORT. 
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FIGURE 30 

REFERRED STATOR TORQUE VS. STATOR PRESSURE RATIO 



T 



T 



T 



T 




78 



FLOW EXIT ANGLE FROM HONEYCOMB VS. REFERRED FLOW RATE 



APPENDIX A 



PREDICTION OF STATOR DISCHARGE ANGLE 
AND STATOR LOSS COEFFICIENT 

1 . Stator Discharge Angle 

The prediction of the stator discharge angle uses methods described 
by Vavra [Ref. 8] and Ainley [Ref. 2]. In these methods, it is assumed 
that the discharge angle is influenced only by blade geometry and Mach 
number, and that flow incidence angles have no effect. 

For values of M < 0.5, a formula given by Vavra [Ref. 8] is used 
for a first approximation, and the effects of blade geometry and Mach 
number are then accounted for by using the method described by Ainley 
[Ref. 2]. A linear variation of a is then assumed for values of M^ 
between M^ = 0.5 and M^ = 1.0. The results of these calculations for 
the blade geometry shown in Fig. 7 are listed in Table III. 

2 . Stator Loss Coefficient By Method of Ainley 

The prediction of stator loss coefficient uses the method described 
by Ainley [Ref. 2]. In this method it is assumed that the loss coef- 
ficient is independent of Mach number and is influenced only by blade 
geometry and flow incidence. The losses are divided into profile 

losses Y which include mixing losses, secondary losses Y , and tip 
P s 

clearance losses Y, _ where 
k 




and 
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Y 

s 



y k = 



A + B 






Values of Y = f(s/C) are given in Fig. 4 of Ref. 2 and values 

P 




are given in Fig. 8. A relation given by Vavra [Ref. 9] is then used 
to /relate the stator loss coefficient to the pressure loss coefficient. 
The results of these calculations are listed in Table IV. 

3. Loss Coefficient By Method of Traupel 

The method of Traupel is one of the most complete methods used 
and takes into account a large number of variables considered to have 
an influence on the loss coefficient. 



£ = £ + C + C + C 

p w r zus 



where £ = Total loss coefficient 

£ = Profile loss coefficient 

P 

= Wall loss coefficient 
= Secondary loss coefficient 
= Damping wire loss coefficient 
The profile loss coefficient is given by: 

C = C 



w 



zus 



'p ’pO * X pO 



• X M * X s + + 5 F 



where £ ^ = Basic loss coefficient 
pO 

• XpQ = Reynolds number correction 



^Traupel, W., n Thermische Turbomaschinen ? Erster Band /* Springer 
Verlag, Berlin/Goettingen, Heidelber, 1958, pg. 269-298. 
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y w = Mach number correction 
M 

X g = Trailing edge thickness correction 

C = Mixing loss correction 
m 

~ Correction for straight blades 
r 

The wall loss coefficient, which is due to the friction loss at the 
annulus of the hub and tip, is given by: 

s. sin 

C w = So ’ x p0 h 

The secondary loss coefficient is given by: 



S “ ^ Z ^rO + C s 

where ip^ = Aspect ratio correction 

C r Q = Basic secondary loss 
= Tip leakage loss 

A summary of the calculations is listed in Table V. 

4 . Loss Coefficient Method By Soderberg 

The method of Soderberg is taken from a paper by Horlock. 
Soderberg bases the loss on blade geometry, aspect ratio, thickness 
ratio, and Reynolds Number, and then relates the loss coefficient to 
the gas deflection angle through the aspect ratio only. He ignores 
such parameters as Mach number and certain non-dimensional gas prop- 
erties . 

The results of a calculation for the converging TTTR nozzles is 
listed in Table VI. 



^Horlock, J. H 
Ins. J. Mech. Sci. 



, "Losses and Efficiencies in Axial Flow Turbines 
Vol . 2, 1960, pg. 59-60. 



w 
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TABLE III 



Prediction of Stator Discharge Angle 



s = 0.8594 in 
a = 0.2050 in 
t = 0.0240 in 



a/s = 0.2385 in 
t/s = 0.0279 in 
For <_ 0.5: 




= 0.9809 

Hence a = cos ^ (0.24318) = 75.92° 

a = a * + 4 ( i) 



where e = 



= 4.5338 
o 



4 (s/e) = 0.76 



Hence a = 76.68 



For M = 1.0 



a (M^ = 1.0) = cos 



-1 



s - t/cos a 



= cos 1 (0.27152) 



For (average) = 0.65 



(average) = 75.90 



[Ref. 8] 



[Ref. 2] 



= 74.25° 



82 



TABLE IV 



Loss Coefficient by Method of Ainley 



= 0 



a 2 = -75.90 



t = 0 

c = 1.4313 
s = 0.8594 
s/c = 0.6004 



Y , 1 = 0.0490 

p(al = 0) 



Y + Y = 
s K 



X + B 



00 



s/C 



2 

cos a„ 
z 

3 

cos a 

m 



tan a + tan a ? 

where tan a = = -1.9906 

m 2 



a = -63.33 
m 

B = 0.5 For non-shrouded wheels 
C L 

— 777 = 2 (tan a- - tan a,,) cos a = 3.5720 
s/C 1 2m 



X = f 



( A 2^ A l) 2 

(1 + I.D./O.D.) 



A 2 A 02 COS a 2 



A 1 A oi cos a l 



itD ^2 h 2 cos a 2 

irD .. h_ cos a, 
ml 1 1 



h l ’ V D ml ' “m2 



Hence X = f 



2 , 2 
cos a^/cos 

+ R hubl /R tip 1 )_ 



= f (0.0321) = 0.0056 



Fig. 4 



Fig. 8 
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K = 0 



Y + Y = 0.0056 (3.5720) 2 cos ( 75.90) = 00471 

cos (-63.33) 

Hency Y = Y + Y +Y, = 0.0961 
p s k 

t/s = 0.0279 

Trailing edge correction factor F = 1.02 
Hence Y = 1.02 (0.0961) = 0.0980 




[Ref. 9] 
page 99 



With (average) 




= 0.7055 (Run 15) 
average 



(, = 0.0756 

s 
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TABLE V 



Loss Coefficient By Method of Traupel 

C n = f(a n a-) = 0.0040 Fig. B3 

pO U , 1 

, - KS, 

+p0 ‘ £ l (Re ' “> 

Re . V . 750 x 1,4313/ 12 . . 5 . 72 x 1() 5 
V 1.564 x 10 ^ 



K = 0 



Hence x o “ 0.92 

p0 

ib = f(Mj = 0.075 for M =0.65 
Ml 1 

= f(l-K) 

1 - K = 1 - ( s \ = 1 - 

cos a^J 






0.8594 



0.0240 
cos 75.9° 



Fig. B5 
Fig. B4 



0.2392 



Hence \b =1.60 
s 

C M = f (1 - K) = 0.0700 
C = f (h/D ) 

r m 



h/D 



m 



0.690 

8.48 



0.0814 



Fig. B6 
Fig. B6 



Hence 4 = 0.0250 Fig. B7 

r 

Thus £ = 0.0040 x 0.92 x 0.075 x 1.60 + 0.0700 + 0.0250 

P 

or ^ = 0.0954 

P 



4 



w 



= 0.0040 x 0.92 x 



0.8594 x sin 75.9 
0.690 



0.0044 



V* B o 



c/h = 



= f (c/h) 

1.4313 = 
0.690 



2.078 
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1.70 



ELg. B9 



Hence = 

and = 1.70 (0.92) = 1.563 




Hence £ n = 0.0350 
rO 

C = f(h/h) = 0 
s s 

Hence c - 1.563 x 0.0350 = 0.0547 
r 



From point 8 Run 15 



Fig. B8 

Not applicable to stator 



4 =0 Not applicable to stator 

zus 

Hence £ = 0.0954 + 0.0044 + 0.0547 = 0.1545 

s 
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TABLE VI 



Loss Coefficient By Method of Soderberg 



5 1 = ( 1 + e*) 0.975 + 0. 



ex _ — 0 
0 



cr^- 75.90 



b = 0.975 
h = 0.690 
t/c = 0.248 

E* = f(Aa) = 0.0700 
Hence 5' = 0.1580 



- 0 - ■ 



Blade axial chord 
Blade height 
Thickness /chord 



Fig. 5 




With (average) = 0.65, 




= 0.7055 (Run 15) 
average 



C = 0.1022 

s 
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APPENDIX B 



PREDICTION OF TORQUE ON CLOSURE 



PLATE AND SHROUD 



1. Torque on Closure Plate 

Because of the proximity of the rotor and its high speed rotation, 
a core of rotating fluid will be set up between the rotor and closure 
plate, and friction will thus transmit some torque to the closure 
plate . 

Schlichting ^ gives a formula for the torque transmitted to a 
stationary flat plate by a rotating parallel plate, assuming the core 
of fluid rotates at one-half the speed of the wheel. 

For turbulent flow: 



This results in a moment of 0.007 ft-lbs at 10,000 RPM and 0.023 ft-lbs 
at 20,000 RPM. 

2 . Torque on Shroud 

A torque will be transmitted to the rotor shroud by the tangen- 
tial component of and by the scrubbing of rotor at the tip. A 

^Schlichting, H., "Boundary Layer Theory," pg. 607 



M = 0.037 pu) 2 R 5 




1/5 



Using values of p and v for air at 520°R and one atmosphere, 



with R = 3.87 in 



-fi 



1 1/5 



M = 0.037 x 0.002378 
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formula given by Schlichting for the torque transmitted between two 
rotating cylinders is: 



2 2 




With = 0 and r = r^ , 

M - T A 2 r 2 - >■ (Sr) (_ 2 * r 2 h ) r 2 




As before, with the core assumed to ratate at one-half the rotor speed, 
and the values for air at standard conditions, the moment transmitted 
is 0.018 ft-lbs at 10,000 RPM and 0.037 ft-lbs at 20,000 RPM. 



8 



Schlichting, 



H. 



5 



"Boundary Layer Theory, " pg. 81 
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APPENDIX C 



PREDICTION OF STATOR DISCHARGE 
PRESSURE DISTRIBUTION 



The equation of motion for absolute flow through a stator, as 
given by Vavra [Ref. 9] is: 

+ VH = V x (V x V) + T Vs + 1 (C.l) 

If steady, frictionless flow is assumed, the equation reduces to: 



VH * V x (V x V) + T Vs (C. 2) 
Each variable in the above equation may then be expressed in an axisy- 
mmetric coordinate system, where 



H \ 3m R 90 1 n 9n 



T Vs = T ( i + ff + i 0 
\ m 3m R 30 n 3n J , 



v 0 

and V x (V x V) = i m -n — 



3(RV Q ) 

9m 



9V 



m 



90 



- V 



n 



9V 



-r^ + V K 

on mm 



+ V 



3V 

— + V K 
3m n n 



+ M v 

R * m 



9V 9 (RV) 
m 0 



30 3m 



+ V 



3V 
n 

30 



3V 

m 




3(RV 0 ) 



9n 



9V 
n 

90 



+ V 1 "T + V K 
m 9n mm 



1 - V + v K 

n [ m j 3m n n_ 



Equating components : 
~ V 



*t . 3H 
m‘ 3m 



0 



R 



9(RV 0 ) 

9m 



9V 



m 



90 



- V 



9V 



m 



a + V K 
9n mm 



+ 



90 



9s_ 

9m 



(C.3) 



V 

n 



9V 

— 2: + V 
dm n 



i 

K 

n 



1 


9H 


V 

m 


9V 

m 


3(RV 0 ) 


V 

+ -JEL 


3V 

n 


3 (RV ) 


R 


96 


R 


30 


9m 


R 


_ 96 


9n _ 



1 9_ 

R 96 



(C.4) 



- . _9H 

1 n * 9n 




8 ( R V 5, 

3n 30 



+ V 

m 



9V 

— H + V 

9n m 




- V 

m 



9V 

— 2. + V 
9m n 





(C.5) 



where m = streamline coordinate 

n = coordinate normal to streamline 

K = curvature of streamline 
m 

K = curvature of normal 
n 

If axisymmetric flow is assumed, all terms with 9/90 vanish and since 
the equation is written along a streamline where = 0 



1 J)H 

R 96 



V 




9(RV 0 ) 

9m 



9(RV 0 ) 

Since V is not zero, — 

m 9m 



0 



(C.6) 



(C.7) 



Then — = — and if adiabatic flow is assumed then — = 0 means H is 
dm dm dm 

constant along a streamline. In addition, - = 0 means the flow is 

dm 

isentropic along a streamline. A consequence of these last two state- 
ments results in isentropic flow or V x (V x V) = 0. The remaining 
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equation is then: 



an v e 9(R V 



3V 



3n R 3n 



+ V 



m + V 2 K + 1^- 



m 9n m m 9n 



(C. 8) 



or 






2 v e 3(R V 3H 

+ V 2 K + — — - + T —■ = 0 

m m R n dn dn 



2 2 
v e + V m 

Employing the relation H = C^T + [ J the following important 



( 



equation results: 

2 



y-(v 2 \+V 2 (2K - -f-] + 2 

dn m ) m\^m C dn J - R 



v a 9(RV 0 ) 



3T„ 



3n 



- 2 C 



p 3n 



\ /, 



+ 2t t-c- 



ds 

3n 



= 0 



(C.9) 



Equation (C.9) is known as the Radial Equilibrium Equation and is a 

linear differential equation with variable coefficients. 

If three additional important assumptions are made, namely that 

there is no temperature or entropy gradient at the stator entrance, 
3(RV q ) 

and that — : = 0 or free vortex flow exists ahead of the stator. 

dn * 

and finally that the streamlines are cylindrical surfaces with K = 0, 

g / 2 \ 

then — { V ) = 0, which means V is not a function of radius. If 
3n V m J m 

— (RV.) = 0 can be maintained in the stator, the same conclusion 
dn o 

can be made for V at the stator discharge, 
m 

The consequence of the above series of assumptions results in a 
set of velocity profiles as shown in Fig. 32. Also imposed on the 
profiles shown in Fig. 32 are the effects of annulus friction and 
secondary flow. 
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(a) Axial 

Component V 

m 



(b) Tangential 
Component V 

D 



(c) Resultant 
Velocity V 



Figure 32 

Predicted Velocity Profiles at Stator Discharge 



The unequal pressure between the convex and concave sides of a 
blade cause trailing vorticies to be formed and are the origin of 
the secondary flow. The vorticies cause a reduction in the tangen- 
tial component near the hub and tip radii. It is evident from the 
above velocity profile in Fig. 32C that the pressure distribution 

must necessarily be a complicated one and not linear. 

9(RV 0 ) 

When the assumption that — = 0 is put in cylindrical coordi- 

□n 

nates, so that: 



B (RV ) 
u 

BR 



B (RV^ tana) 
BR 



= 0 



or 

V R d . (tana) . + v tana _ 0 
a 9R a 

hence 

8 (tana) 9a tana 

SR 9R R 



(C.10) 
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But for 



= 0 the blades must themselves be radial, and their 



8 (RV ) 

“3R 

geometry is given by Vavra [Ref. 9] as 



3a _ tana 
3R R 



(C.ll) 



Since equations (C.10) and (C.ll) are in opposition to each other, 

the blades cannot everywhere be radial, and thus the blades will 

exert a force on the flow in a radial direction, which will produce 

curved stream surfaces. As a result, the assumption of axisymmetric 

3 (RV ) 

stream surfaces is in error, and in addition — - = 0 cannot be 

d R 

maintained in the stator, so free vortex flow is not possible. Thus 
2 



3V 



V 3(RV ) . .. , - , . r . . . - 

a « __u u and this can be solved if it is assumed there 

3R 1 R 3R 

is constant total enthalpy and flow angle ahead of the stator. The 

solution is given by Vavra [Ref. 8] which is: 



. 2 

sm a, 



V = V ^ 
a am 




where V = axial velocity at the mean raduis . Thus V = f(R) and 
am y a 

the resulting velocity or pressure distribution becomes extremely 
complicated. 
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Printout of Programs FLOCAL and TTTR for use of IBM 360/67 
computor . 
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WRITE(6,20)PAMB,N, GAM 

WR I TE ( 6,21 ) 

DO 13 1 = 1, N 

WRITE (6, 22 ) I , WFLI I ) , QKN ( I) , RE ( I ) 
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